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Abstract 

The Makeenko-Migdal loop equation is non-linear and first order in the area deriva- 
tive, but we show that for simple loops in QCD2 it is possible to reformulate this 
equation as a linear equation with second order derivatives. This equation is a bound 
state Schrodinger equation with a three dimensional Coulomb potential. Thus, loop 
dynamics leads to a surprising new picture of confinement, where this phenomenon 
is due to a (bound state) localization in loop space, with the Wilson loops decaying 
exponentially outside a characteristic radius. 
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The Makeenko-Migdal loop equation [1] is an exact consequence of the non-hnear dy- 
namics in A'^ = oo QCD, and is therefore also non- linear. It is of first order in the area 
derivative. However, although the equation is closed, loop space is enormous, so it is not 
easy to gain insight into how dynamics look like in this space. In the present letter we 
restrict ourselves to simple loops with an arbitrary number of windings n corresponding to n 
quarks propagating along the boundary of the loops. We then reformulate the loop equation 
for these simple loops in QCD2 starting from the Makeenko-Migdal equation. The result is 
an equation which is linear but of second order in the derivatives. This is obtained from 
the Hopf equation in QCD2, which can be generalized to a functional equation in higher 
dimensions, so the two dimensional result may be generalizable to higher dimensions. In any 
case, our linear equation represents a new picture of confinement as a bound state in a three 
dimensional space. The interaction between the n quarks is given by the three dimensional 
Coulomb potential in loop space! This shows that the dynamics in loop space is rather 
surprising and is not in a simple correspondence with perturbative expansions . 

For simple loops with a number n of windings it is possible to derive an equation for the 
spectral density of eigenvalues 



Pc{9) 



27r 



1 + 25^ W^''\C)cos{ne) 



n=l 



(1) 



where W^^\C) is the Wilson loop for the simple curve C winding n times. The equation is 
[2] 
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where P means the principal value of the integral. If we introduce the resolvent (i.e. the 
Hilbert transform) 
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we can rewrite eq. (2) as 
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This equation is non-linear, as expected. 

In two dimensions eq. (4) reduces to the partial differential equation [2], 
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de 
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where it was used that in D = 2 the Wilson loops depend only on the area A, which is 
measured in units of QqN . Eq. (5) is called the Hopf (or Burgers^) equation. In [3] it is 
shown that this type of equation, or a slight generalization, occurs for any multiplicatively 



^Although the viscosity term usualy occuring in Burgers eq. is absent here. 
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free family of random variables [4] . Eq. (4) is valid in all dimensions and is a Hopf equation 
in loop space. 

We shall now show that for the boundary conditions 1^^")(0) = 1 the non-linear first 
order Hopf equation (5) reduces to a linear second order equation. Let us use the notation 
/(O, 9) — fo{9). By means of the method of characteristics eq. (5) can be reformulated to 

f{A,e) = MaA,e)), aAO)^d-AMaA,e)), (6) 

and the various derivatives can be expressed as 
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Here the prime means differentiation with respect to ^. From (7) we obtain 
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and 
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By comparison of eqs. (8) and (9) we see that we can form a linear second order differential 
equation of the form 

with a and b to be determined, provided 

a/o (0 = -2/o(0/o(0, (11) 

and 

To solve, insert eq. (11) in (12), which gives 

fo{^{A, 9)) = f{A, 9)^Vb coth(V^ ^A: 9) /a). (13) 

By inserting this result we see that eqs. (11) and (12) are satisfied without any restrictions 
on a and h. 

It now remains to see whether the result (13) fits with QCD2, where 

PA=o{9)= f; 5(^-27rn), (14) 
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SO from eq. (3) we get 
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If we compare this with the solution (13) we see that QCD2 satisfies the hnear second order 
differential equation (10) provided a = 1 and Vb = i/2. Hence the differential equation (10) 
becomes 



This linear equation can thus replace the non-linear Hopf equation (5) when the QCD2 initial 
value (15) is used. 

It should be mentioned that in a similar manner one can produce other second order 
linear versions of the Hopf equation. For example, instead of eq. (16) we could study 

^?^.(.+,.)^. (u) 

which leads to a cubic equation for the initial value function 

c/o(0 - 1/3 MO' = di + const. (18) 

This equation has no relation with QCD2, and the linear differential equation (17) must 
therefore be rejected. Only eqs. (10) and (16) are relevant. 

As shown in ref. [3] the Hopf equation (or simple generalizations) occurs in any mul- 
tiphcatively free family of random variables. Thus the method used above to rewrite the 
Hopf equation in terms of equations which incorporate the initial conditions can be used to 
obtain a linear replacement of the relevant Hopf equation. The point is then to chose the 
linear equation such that the initial function is physically relevant, which may not always 
be possible. 

Let us now return to the linear QCD2 equation (16). Suppose this is all we know, then 
how do we solve it? The simplest way is to use the factorization f{A, 9) — fi{A)f2{9). Eq. 
(16) can be separated in the two variables, and the resulting equation for f2{d) is trivial 
and has the solution f2{d) = exp(A^), where is the separation constant. Since f{A,6) is 
defined to be periodic in 9 with period 2n, it follows that A should be an integer times i. 
For each such A the function /i should be proportional to a W^^\A). Referring to of eqs. 
(1) and (3) the separation method gives 

f{A,9) + £ , (19) 

and the Wilson loops should satisfy the separation equation 



dA'^ A dA 4 \A 
This equation has the solution 

W^''\A) = -L^_l(nA)e-"^/^ (21) 
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where the L's are the Laguerre polynomials of type 1, and where we used the boundary 
condition W^'^")(0) = 1 to fix an arbitrary constant. Eq. (21) is the well known solution for 
winding loops [5] in QCD2. 
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Eq. (20) can of course be derived directly from the known solution (21). The reason we 
have started from the Hopf equation (5) is that it can be generalized to higher dimensions 
as is seen in eq. (4). In D = 3 or 4 it is not true that the resolvent fc{0) only depends 
on the area, even for simple loops. However, one may speculate that area dependence sets 
in asymptotically for sufficiently large loops. On the other hand, from the works in refs. 
[6] it follows that the loop equation in D = 4 has features quite different from its D = 2 
counterpart. It would be an interesting project for future research to try to invent a general- 
ization of the method of characteristics (leading to eqs. (6) and (7) in two dimensions) to the 
functional version of the Hopf equation valid in higher than two dimensions, incorporating 
the complexities discussed in refs. [6]. 

In a recent work by Narayanan and Neuberger [7] they have considered the spectral 
density in four dimensions and found that it is similar to the density in two dimensions. This 
approach refers to Wilson loops in which the Coulomb-like contribution has been smeared in 
order to avoid divergences, thereby making the D — A case similar to the D — 2 case. If the 
the functional Hopf equation (4) could be similarly smeared it would probably be similar to 
its two dimensional version for large smeared loops. 

Eq. (20) is the (zero angular momentum) radial Schrodingcr equation for the hydrogen 
atom if we identify the radial variable with the loop quantity A and the charge with n. 
Since n quarks are running around on the boundary, this is a quite natural identification. 
Similarly the "bound state" energy is — n^/4. The point A = 4 is essentialy the "Bohr 
radius" of the atom. Inside this radius the "wave function" W^'^^ oscillates as a function of 
A, whereas outside this quantity is exponentially damped. Although confinement is often 
considered trivial in two dimensions, the "wave functions" II/(") behave in a rather non- 
trivial manner, especially for larger n's. Similar behavior does not occur in the Abelean case 
in two dimensions. Since there exist the hope that the Hopf equation can be treated by a 
functional version of the method of characteristics in three or four dimensions, the bound 
state picture of confinement may be valid also in higher dimensions. 

Loop space dynamics thus results in a rather surprising picture of confinement, where 
the area acts as a radial variable in a three dimensional space in the sense that the force 
between the quarks with charges n is governed by the three dimensional Coulomb potential. 
The winding Wilson loops W^^^\A) act as wave functions satisfying a Schrodinger equation. 
Inside the localization radius A = 4 these wave functions oscillate (more and more when n 
increases) and outside A = 4 they decay exponentially. These three dimensions are genuinely 
a product of the loop dynamics, which seems to have forgotten about the two dimensions 
where it lives. It would be interesting to investigate other more complicated loops to see to 
what extent this picture of confinement is universal. 

One can speculate that the reason for the localized bound state picture is that QCD2 
basically is a family of multiplicatively free random variables [3] , which is in analogy with 
Anderson localization of the wave function in the presence of random fields in quantum 
mechanics. It also could well be that a similar phenomenon occurs in higher dimensions due 
to the expected random fields in the QCD vacuum. 
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